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Abstract 

We briefly review the five possible real polynomial solutions of hypergeometric differential equa- 
tions. Three of them are the well known classical orthogonal polynomials, but the other two are 
different with respect to their orthogonality properties. We then focus on the family of polynomials 
which exhibits a finite orthogonality. This family, to be referred to as the Romanovski polynomi- 
als, is required in exact solutions of several physics problems ranging from quantum mechanics and 
quark physics to random matrix theory. It appears timely to draw attention to it by the present 
study. Our survey also includes several new observations on the orthogonality properties of the 
Romanovski polynomials and new developments from their Rodrigues formula. 
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I. INTRODUCTION 



Several physics problems ranging from ordinary-and supersymmetric quantum mechan- 
ics to applications of random matrix theory in nuclear and condensed matter physics are 
ordinarily resolved in terms of Jacobi polynomials of purely imaginary arguments and param- 
eters that are complex conjugate to each other. Depending on whether the degree n of these 
polynomials is even or odd, they appear either genuinely real or purely imaginary. The fact 
is that all the above problems are naturally resolved in terms of manifestly real orthogonal 
polynomials. These real polynomials happen to be related to the above Jacobi polynomials 
by the purely imaginary phase factor, i", much like the phase relationship between the hy- 
perbolic and the trigonometric functions, i.e. sin ix = i sinh x. These polynomials have first 
been reported by Sir Edward John Routh [1| in 1884, and then were rediscovered within the 
context of probability distributions by Vsevolod Romanovski 2f| in 1929. They are known 
in the mathematics literature under the name of "Romanovski" polynomials. 

Romanovski polynomials may be derived as the polynomial solutions of the ODE 

(l + x^)^ + t(x)^ + Ai?(x) = 0, (1) 
with t{x) a polynomial, at most a linear, which is a particular subclass of the hypergeometric 



differential equations 



^ . Other subclasses give rise to the well known classical orthogonal 



polynomials of Hermite, Laguerre and Jacobi jj], j^. Romanovski polynomials are not so 
widespread as the others in applications. But in recent years several problems have been 
solved in terms of this family of polynomials (Schrodinger equation with the hyperbolic 
Scarf and the trigonometric Rosen- Morse potentials la, LDl, Klein-Gordon equation with equal 

fi h 

vector and scalar potentials [3], certain classes of non-central potential problems as well [9|) 
and so they deserve a closer look and be placed on equal footing with the classical orthogonal 
polynomials. 

In this context, our goal is threefold. First of all, it is to establish the orthogonality 
properties of these polynomials. This is achieved by the same methods as for any other 
hypergeometric differential equation. Our second goal is to explain their use as orthogonal 



eigenf unctions of some Hamiltonian operators. Third, Eq. ([T]) has been described in 10|] as 
a complexification of the Jacobi ODE, a general expression that can be written as 



2.d^P{x) , ,,^^dP{x) 



;i - x')—^ + ^(^)^ + ^n^) = 0, (2) 



where t{x) is again an arbitrary polynomial of at most first degree, but not necessarily the 
same as in Eq. ([T]). If that were the case, solutions to Eq. ([T]) would be the complexification 
of the solutions to Eq. ([2]), that is, the complexification of the Jacobi polynomials. Hence, 
our final goal is to clarify this relationship. 

We deal with all these issues in the following way: In Section [IT] we give a classification 
of hypergeometric equations placing Eq. ([T]) among them. Next, in Section IIIII we show 
some expected properties of the functions as solutions of a hypergeometric ODE 

such as: being indeed polynomials, recurrence relations; and the absence of another, namely 
general orthogonality. In Section [IV] we compare the polynomials Ri^''^\x) with the com- 
plexified Jacobi polynomials. In Section |V] we show some examples of physical problems 
whose solutions lead to Romanovski polynomials. Section IVII sheds light on some peculiar- 
ities of orthogonal polynomials as part of quantum mechanics wave functions. In the final 
Section fVIII we summarize our conclusions. 

II. CLASSIFICATION OF HYPERGEOMETRIC DIFFERENTIAL EQUATIONS 

A hypergeometric equation [s^ is an ODE of the form 

s{x)F"{x) + t{x)F\x) + \F{x) = 0, (3) 

where the unknown is a real function of real variable F : U H, where W C R is some 
open subset of the real line, and A G R a corresponding eigenvalue, and where the functions 
s and t are real polynomials of at most second order and first order, respectively. Here the 
prime stands for differentiation with respect to the variable x. This class of ODEs is very 
well known both from the mathematical and the physical points of view. From the mathe- 
matical one, many properties that its solutions exhibit make them interesting in their own 

Ann 

right. For instance, the classical orthogonal polynomials [3[, [12j (Hermite, Laguerre 



and Jacobi polynomials, the latter including as particular cases Legendre, Chebyshev and 
Gegenbauer polynomials) are solutions of particular subfamilies of hypergeometric ODEs. 
From the physical point of view, many of the exact solutions to the eigenvalue equation of 
a quantum mechanical Hamilton operator lead to an equation of the hypergeometric kind: 
harmonic oscillator. Coulomb potential, the trigonometric Rosen-Morse and Scarf potentials, 
hyperbolic Rosen-Morse and hyperbolic Scarf potentials. 
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As to our goal, the mathematical properties we are interested in are the following (refer 
to for a detailed study and proofs of these statements). The leading property, which gives 
the differential equation its name "hypergeometric," is that if F{x) is a solution to Eq. ([3]), 
then the derivative F'{x) is a solution to another hypergeometric equation that is closely 
related to the former: 

s{x){F'{x))" + t^^\x){F'{x)y + AWf'(x) = 0, (4) 

where t^^\x) = t(x) + s'{x) and A*-"*^^ = A + t'{x). Iteratively, it is easy to show that the mth 
derivative, F^"^^ (x) is a solution of 

+ + A('")f('")(x) = 0, (5) 

where now t^'^^x) = t(x) +ms'{x) and A*^™-* = X + mt'(x) + ^m{m — l)s"{x). The next result 
is that, for any n G {0, 1, 2, . . . }, there exists a polynomial Fn{x) of degree n, together with 
a constant A„ which satisfy Eq. ([3]). The constant is given by 

Xn = -n(^t'{x) + ^{n-l)s"{x)^ . (6) 

The last result, together with the former, tells us that -F„(x) and its derivatives, Fn ^ (3;) , are 
solutions to similar equations. By means of a weight function, it is possible to write down a 
formula which gives all these polynomials at once. A weight function w{x) associated with 
Eq. ([3]) is a solution of Pearson's differential equation 

[s{x)w{x)]' = t{x)w{x), (7) 

that assures the self-adjointness of the differential operator of the hypergeometric ODE. 
Then, the generalized Rodrigues formula gives the mth derivative of the polynomial Fn{x) 
as 

1 Jn—m 

= Ar„„— -— - — [w(x)s(xr] , 

^ ' w{x)s{x)-^ dx'^-'^^ V ; V ; J, 

< m < n, 

where Nnm is a normalization constant. This constant is related to the coefficient a„ of the 
term of degree n in the polynomial Fn{x) by the expression 



which is vahd for < m < n — 1 and n > 1. Equation ([H]), with m = 0, gives the classical 
Rodrigues formula 

1 d"- 

= iVn^y^ [wix)sixr] , (10) 
w[x) dx"- 

where we have identified Fn^\x) = Fn{x) and A^„o = Nn- 

If the functions s{x) and w{x) satisfy yet another condition, namely both being positive 
within an interval (a, b) and 

lim s{x)w{x)x'' — lim s{x)w{x)x'' = , (11) 

X — >a X — >b 

for any nonnegative integer /, then the family of polynomials is orthogonal with respect to 
the weight function w, i.e. 

rb 



/ w{x)F^{x)Fn{x) dx = {fnf^mn, 
J a 



(12) 

Vm,n G {0,1,2,...}, 

where /„ is the norm of the polynomials. Hence, all hypergeometric ODEs admit a family 
of polynomial solutions. But this family is not orthogonal for all hypergeometric ODEs. 

The fact that a solution F{x) and its derivatives F^'^\x) obey hypergeometric ODEs 
with the same coefficient s(a;), Eqs. ^ and ([S]), suggests a classification in terms of the 
polynomial s{x). Moreover, a classification according to the roots of s{x) has proved useful 
and provides a characterization of the solutions j3|,jj],j^. There are five classes in this 
scheme, as s{x) may be a constant, a first degree polynomial or a second order one with 
two distinct real roots, one real root or, finally, two complex conjugate, not real, roots. In 
addition, it is useful to note that an affine change of variable (i.e., x— >ax + 6, 07^ 0) 
does preserve the hypergeometric character of Eq. ([3]) and the kind of roots of polynomial 
s{x). Then, in each class, we may consider only a canonical form of the equation, to which 
any other can be reduced by an affine change of the independent variable. 

1. Polynomial s{x) is a constant: 

We take as canonical form 

H"{x) - 2axH'{x) + XH{x) = 0, (13) 

where a G R is an arbitrary constant, i.e., we have here a one-parameter family of ODEs. 
We call it generalized Hermite equation (the equation with a = 1 is called Hermite equa- 
tion). The polynomials are a generalization of Hermite polynomials, denoted {ifn""*}. 



n E {0, 1, 2, ... }. The weight function is 

= e""^'. (14) 



For a > the additional conditions for orthogonahty, Eq. (fTTI) . are fulfilled in the interval 
(— oo, oo), hence we get an orthogonality relation: 

(15) 

Vm,n e {0,1,2,...}, a > 0. 

2. Polynomial s{x) is of the first degree: 
The canonical form of the ODE is 

xL"{x) + t{x)L\x) + XL{x) = 0, (16) 

which we call generalized Laguerre equation. The first-order polynomial t is still arbitrary, 
so this is a two-parameter family of ODEs. If t{x) is written as t{x) = —ax + f3 + 1, with 
a,/3 G R, the parameters (actually, Eq. ( |T6|) is called associated Laguerre equation in the 
case a = 1, and Laguerre equation if a = 1 and (3 = 0), then the weight function is 

w{x) = a;^e-"^ (17) 

and the polynomials are written {lI^'^^}, n G {0,1,2,...}. li a,f3 > 0, the condition of 
Eq. (|TT1) is fulfilled and one gets orthogonality in the interval [0, oo) as 



(18) 



■^^^ ■^rn'^\-^)'^n dx — (In) 5mn; 

Vm,n G {0,1,2,...}, a,/? > 0. 

3. Polynomial s{x) is of the second degree, with two different real roots: 
The canonical form of the ODE is 

(1 - x^)P'\x) + t{x)P\x) + \P{x) = 0, (19) 

which is known as Jacobi equation. It is customary to write the arbitrary polynomial t{x) 
in the form t{x) = f3 — a — {a + (3 + 2)x, where a, /5 G R are the parameters. Then, for each 
pair [a, (3), the Rodrigues formula defines a family of polynomials, the Jacobi polynomials, 
denoted {Pt'^^}, n G {0, 1, 2, . . . } with the weight function given by 

w{x) = {l-x)''{l + x)^. (20) 



If parameters a and (5 satisfy > —1, the additional condition of Eq. flTTl) is fulfilled in 
the interval (—1, 1), so there is an ort honor malizat ion relation: 

' (1 - xYil + x)^Pi"''3)(x)P^'3)(x) dx = (J9„)25^„, 

(21) 

Vm, n G {0, 1, 2, . . . }, a, /3 > -1. 

Some particular cases received special names: Gegenbauer polynomials if a = /3, Chebyshev 
I and II if a = /3 = ±1/2, Legendre polynomials if a = /3 = 0. 

4. Polynomial s[x) is of the second degree, with one double real root: 
We choose as canonical form of the ODE 

x'^B"{x) + t{x)B'{x) + XB{x) = 0, (22) 

If the arbitrary first order polynomial is written as t{x) = {a + 2)x + (3, with a, /3 G R the 
parameters, then the weight function is 

=a;"e-i. (23) 

We write the polynomials as n G {0, 1, 2, . . . }, which are called Bessel polynomials 

[isl (they were also given under type V in Ref. and classified in Ref. j4|). There is 
no combination of any particular values of the parameters and any interval which satisfies 
Eq. (ITT!) , so neither of these families is orthogonal with respect to the weight function (l23ll . 

5. Polynomial s{x) is of the second degree, with two complex roots: 
The canonical form of the ODE for this case is chosen as 

(1 + x'^)R"{x) + t{x)R'{x) + XR{x) = 0, (24) 

which is the one studied in and Q [isl. A note of caution: the solutions introduced 

in 6] and [3] seem to be different, but this is due just to a different form given to the 
arbitrary polynomial t{x). A careful review shows that both papers are dealing with the 
same ODE, namely Eq. (IMl) . so the solutions must be the same up to a constant factor. 
Writing the polynomial t{x) as t{x) = 2(3x + a, with a, (3 G R, we have again a two- 
parameter family of ODEs with their respective families of polynomials which we denote 
{rI^''^^}, n G {0, 1,2, . . . }. With this notation (which differs slightly from |6|,Q,[l6]) the 
weight function is 

w^''^^\x) = (1 + x')^-V"™*"'". (25) 



Upon comparison with Romanovski's original work [2.], we conclude {rI^'^^} are the Ro- 
manovski polynomials. In Section IIIII we study the properties of these polynomials. 
In Ref. [4 1 polynomial solutions of linear homogeneous 2nd-order ODEs 

s{x)y"{x) + t{x)y'^{x) = n[{n - l)e + 2e]?/„(x), 

s{x) = ex^ + 2fx + g, t{x) = 2e + 7, (26) 
n G {0,1,...}, e,/,^,e,7 G R, 
are classified upon substituting the finite power series 



ynix)=J2^ix + c)\ 



k\ 

k=0 



(27) 



an,n 7^ 0, C G C , 

and analyzing the resulting recursions among the coefficients. No other solutions other than 
the polynomials given above are found. Their orthogonality properties are derived by means 
of the spectral theorem (of Favard 

To sum up, all hypergeometric equations fall into one of these five classes. Three of them 
give rise to the very well studied classical orthogonal polynomials (Jacobi, Laguerre and 
Hermite) or a slight generalization of them. A fourth one has not attracted much attention 
due to, we guess, the lack of general orthogonality. Finally, a fifth class is the family of 
ODEs we are dealing with here. 



III. DEFINITION AND PROPERTIES OF ROMANOVSKI POLYNOMIALS 

We focus now on the Romanovski polynomials rI^'^^ and study some well known, and 
some new properties they have. We start by writing down the explicit Rodrigues formula, 
Eq. ffTOl) . for this case, that we take as their definition. For each G R and each 
n G N = {0, 1, 2, . . . } we define the function Rih'^\ by the Rodrigues formula 

1 d" 

where 

w(°'^)(x) = (l + x^)^-^"™*"'", (29) 
is the weight function, same as in equation (l25i) . and 

s{x) = 1 + x^ (30) 



is the coefficient of the second derivative of the hypergeometric differential equation 
Notice that we have chosen the normahzation constants A^^ = 1, which is equivalent to 
make a choice of the coefficient of highest degree in the polynomial, as given by equation 
([9]), which takes the form 

n-l 



n = ^ IT Win -k)+ n{n - 1) - k{k - 1)] 



^= .Vo ' ' (31) 

n>l. 

Notice that the coefficient a„ does not depend on the parameter a, but only on (3 and, 
for particular values of /3, a„ is zero (i.e., for all the values {3 = ^^^~2(^„"l"fc"^''''' where k = 
0, . . . ,n — 1). This observation poses a problem that we will address somewhere else. For 
later reference, we write explicitly the polynomials of degree 0, 1 and 2 

4"''^(x) = 1, (32) 
i?S"'^)(x) = (w'^''^^\x)s{x) + /(a;)«;("'^)(x)) = &^^\x) = 2(3x + a, (33) 

Rt^\x) = , , j-[g'(a:)w("-^)(x) + 2g(x)gYx)w("'^)(x)] 

= {s{x)w^'''^\x){&^^\x) + s\x))) 

w(°'^)(x) dx^ ' 

= {2x + &^^\x))&^^\x) + (2 + t'("''^)(a;))s(x) 

= (2/3 + l)(2/3 + 2)a;^ + 2(2/3 + l)aa; + (2/5 + a2 + 2) (34) 

that derive from the Rodrigues formula ( l28l) in conjunction with Pearson's ODE (I7l). 

The whole set of Romanovski polynomials is spanned by the three parameters a, /3 and 

n 

{RI^^^^ : a,/3 G R, n e N}. (35) 

In order to study their properties we have found it useful to classify them in families. The 
properties are stated for each family. We have found two different classifications in families 
of different kinds which share different properties, so we distinguish them in the following 
two subsections. 
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A. The 7^("'^) families 

The family T?.'-"'^-' contains the polynomials with fixed parameters a and p. 

^ia,p) ^ |^(a,/3) . ^ £ N}. (36) 

This family has one polynomial, and only one, of each degree; two different families do 
not share any polynomial in common and the union of all of them gives the whole set of 
Romanovski polynomials (i.e., they form a partition of this set). 

The first property of one of these families is that which led to their construction: the 
family TZ^°''^\x) comprises all the polynomial solutions of the hypergeometric differential 
equation 

(1 + x'^)R"{x) + {2(3x + a)R'{x) + XR{x) = 0, (37) 

where A is a constant which, for the solution E!h'^\x), is given by A„ = —n{2(3 + n — 1). 

Other characteristic properties of classical polynomials are present in the 71^°^'/^^ family 
too, such as a differential recursion relation and an expression for a generating function in 
closed form. 

The differential recursion relation is obtained from the Rodrigues formula, Eq. (l28l) . for 
the polynomial Rn+i{x) (superscripts (a,/3) omitted for clarity). 

1 fjn+l 

= ^d?^ H^H^r^] . (38) 

Then, because of the very definition of the weight function, Eq. ([7]), it is easy to see that 

[w{x)s{x)''+^]' = w{x)s{x)''[2{P + n)x + a]. (39) 

Upon substitution in Eq. (155]) and a straightforward derivation one gets 

Rn^iix) = ^ ([2(/5 + n)x + a]^ [w{x)s{xr] + 2n{(3 + n)^ [w{x)s{xr]^ . (40) 

In the first term, the formula for Rn{x) is recognized, while in the second term its derivative, 
-R^(x), appears (by means of Eq. ([8]) applied to the present case). The result, which makes 
use of Eqs. (|9]) and fl3T]) . is the following differential recursion relation 

2(/5 + n)(l + a:') "^^ = (2/3 + n - 1) (R!^:f^{x) - [2{(3 + n)x + a]R^;:^^\x)) . (41) 
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An integral representation of the Romanovski polynomials is obtained by means of the 
Cauchy's integral formula. As the weight function can be extended to the complex plane, 
where it is analytic except at points ±i, we can use Cauchy's integral formula to get 

^-L/-'°'^^'M-(--)"d,. (42) 

where x is real but z is a complex variable, and 7 is a closed curve in the complex plane, 
enclosing point x, but not ±i. Substituting this equation into the definition of Romanovski 
polynomials, Eq. fl28|) . we get 

2mw('-^^){x)dx- z-x ^ ^ 

The n derivatives with respect to x can be easily performed under the integral sign, giving 
rise to the following integral representation: 

^""^ ~ 27nw(-'^){x) {z-xr+^ ^ ' 

n 

This representation is useful in calculating the generating function, as explained in [3]. A 
generating function, i?*^"'^) (x, y) , of the family 7^(°'^) is a function that is analytic in the 
variable y and whose Taylor expansion in the variable y has the form 

00 1^ 

R^-^^\x,y) = Y^y-Rt'\x). (45) 

A:=0 

Upon substitution of Eq. (jH]) in previous equation, and interchanging the order of integral 
and sum signs (which is allowed since the function is analytic), we get 

' k=0 



The sum is a geometric series, which can be summed as 

Y-(yMy = I = (47) 

"^yz-xj I-MM. z-x-ys{zy ^ ' 



k=0 



z—x 



provided |^^| < 1. In this case, the expression for the generating function results in an 
integral which can be easily evaluated by the method of residues. 

Ri^,^)u y) - I f ^^""^(^) d. = ' Res ( "^"''^("^ z,) (48) 

^ '^^ 2mw^'^'P){x) J.z-x-ys{z) w^^^P){x) \z - x - ys{zy ^ ) ' ^ ' 
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where zi is one of the roots (the one closer to x so it is the only one enclosed by 7) of the 
second order polynomial in z in the denominator: z — x — ys{z) = —yz'^ + z — {x + y). The 
residue is 

\z-x-ys[z) J -4?/(x + ?/) 

where zi = -^{1 — a/I — 4?/(x + y)). Direct substitution gives the final form for the gener- 
ating function: 

^ (2?/A- 1) (1 + x2)/3-ie-"^°t"'^ ' ^ ^ 

where 



^ l- Vl-4y(x + ,) 

2!/ ^ ' 

The next issue to study is the orthogonality properties inside one family. In contrast with 

the classical orthogonal polynomials, these families are not orthogonal with respect to the 
weight function w^°'''^\x) in the natural interval (—00,00) as the following counterexample 
shows. 

Let us consider the family of polynomials with a = and (3 = and, within it, the 
polynomials of degree and 2 which, upon substitution in Eqs. (l32l) and (l34l) . result 

i?r)(x) = 1, (52) 
40'0)(a;) = 2(x2 + l). (53) 

Then, an integral of the form of that in Eq. (JT2l) . in the interval (—00, 00), takes the form 

/oo roo 
w(°'°)(a;)4°'°^(x)4°'°^(x)dx = 2 / dx, (54) 
-00 J —00 



which does not converge. Notice that this poses also a problem with the normalization: 
polynomial R'2'^\x\ for instance, is not normalizable. The most that can be said is the 
following theorem jj],!^! that we state and prove, which establishes (the so-called finite) 
orthogonality among a few of the polynomials in the family. 



If Rm'^\x) and rI^''^\x), m ^ n, are Romanovski polynomials of degree m and 
n respectively, then 

/oo 
ix)Rl^^^^ (x)i?i°'^) (x) dx = (55) 
-00 

if, and only if , m + n < 1 — 2/3. 
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The proof follows the steps of the usual proof of orthogonality of the polynomial solutions of 
hypergeometric equations. Express Eq. ( |3711 in the adjoint form with the aid of the weight 
function (it is possible thanks to the Pearson's relation in Eq. ([7])) for Rm{x) (in the proof 
we drop the superscripts {a, (3) as they are fixed): 

[w{x)s{x)R'^{x)]' + Xmw{x)Rm{x) = 0. (56) 

Multiply by Rn{x), and do the same as above but interchange m and n. Then subtract and 
integrate, which leads to 



{Rn{x)[w{x)s{x)R'^{x)]' — Rm{x)[w{x)s{x)R'j^{x)]'} dx 

/oo 
w{x)Rm{x)Rn{x) dx = 0. (57) 
-oo 

The second integral is the orthogonality integral of Eq. f|T2|) . The first integral, upon in- 
tegration by parts, vanishes except for the boundary terms, so we are left with (because 

Xm 7^ Xn) 

poo 2 

/ w{x)Rm{x)Rn{x) dx = _ {w{x)s{x)[Rm{x)K{x) - R'^i^) Rnix)]}"^^ . (58) 

J —oo '^m 

The boundary term must be evaluated as 

lim M(x) - lim M{x), (59) 

X— »oo x— >— oo 

where 

Mix) = ((l + a;2)V"^°*"^[/2„(a;)i?;(x)-i?:„(x)i?„(x)]). (60) 

Both limits must exist separately. Notice that the limit of the exponential factor is a positive 
constant (1 or e~"'^). The other factors behave as a power of x with exponent (2/5+m+n — 1). 
Then, both limits exist if, and only if, this exponent is negative, in which case both are zero, 
as claimed. 

With a similar argument the following is proved, but we omit the explicit proof: 

For the family of polynomials 7?.*-°''^^ only the polynomials R^'^'^ [x] with n < —(3 
are normalizable, i.e., the integral 

w'^'''^\x){R^^^''\x)f dx, (61) 



converges. 
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So, for the family 7^("'^) only the polynomials in a finite subset are normalizable, and only 
a finite number of couples are orthogonal. This type of orthogonality is sometimes called 
finite orthogonality. 



B. The Q(°'^) families 

In order to calculate explicitly a polynomial i?!"'^'' (x) one has to evaluate the nth deriva- 
tive £^ (w(°'^)(x)s(a;)"), step by step: first £ (u;(°''3)(a;)s(x)") , then ^ (M;("'^)(a;)s(x)") , 
and so on. It is useful to notice the following relation holds 

= (62) 

as can be seen directly in the definitions of w^°''^\x) and s{x), Eq. (!25l) . Thus, it is enough 
to have an expression for the z/th derivative of In looking for it one soon discovers 

the structure of these derivatives: the z/th derivative can be factorized as w{x)s{x)~'^ and 
a polynomial of degree u. For instance, it is easy to calculate directly the first derivative 
(compare with Pearson's ODE ([7])) 

r«;("''^)(a;)l = w^"'^\x)s{x)-^\2{f3 - l)x + a]. (63) 

dx 

By induction, it is easy to prove the statement in general: for fixed a and f3, the z/th 
derivative, z/ G N, of has the form 

d^ 



dx'' 



[w("'^)(x)] = g^(x)w("'^)(x)s(x)-^ (64) 



where Qu{x) is a polynomial of degree z/. The case z/ = is trivial, and the case z/ = 1 is 

shown in Eq. fl63l) . The step from case z/ to z/ + 1 is a straightforward derivation. So, for 
given a, (3 and v we define the polynomial qI"''^'' as 

- . (65) 

Using Pearson's ODE ([7]) as above the first three polynomials for fixed a and {3 are 

Q^^)(x) = l, (66) 

Q(°''5)(a;) = 2(/3-l)x + a, (67) 

Q^^'^\x) = 2{f3 - l){2f3 - 3)x2 + 2a{2f3 - 3)x + 2{f3 - 1) + a^. (68) 

14 



By looking at these expressions, and upon comparison with the first three R^'^^ {x) polyno- 
mials, one readily suspects some relation between the polynomials Qv{x) and Rv{x). This 
relation indeed exists and is the following. 

The polynomial Q^^'^\x) is a Romanovski polynomial, specifically 

Q^^) = R^;;^^-''\ (69) 

The proof is a straightforward manipulation of Eq. (l28l) . applied to '^^ with the aid of 

Eq. ([62]): 



" ^ ' w("./3-'')(x) dx'^ ^ V / V / J 



1 d'^ 



(70) 



w(°'^)(x)s(x)-'" dx' 

But the last term is precisely the definition of Qif'^^x), Eq. (l65l) . 
We now define the family as 

Q(°'^)(x) = {g(,°'^) : 1/ G N} = {i?^,""^-'^) : z/ G N}. (71) 

This family contains Romanovski polynomials with one fixed superscript and the other 
running with the degree. As in the case of the IZ families, the family contains one, 

and only one, polynomial of each degree; two different families do not share any polynomial 
in common and the union of all the families gives the whole set of Romanovski polynomials 
(i.e. they constitute another partition of this set). 

We now study the properties of the polynomials in the family Q("'^). 

In the first place we give some Rodrigues-looking expressions for the Q[f polynomials. 
The polynomials in the family obey the following formulas: 

Q'^''\-) = [w''"'-'\x)sixr] , (72) 

Qi^'^'i-) = [w^''^'\^)s{x)-^Qr\x)] . (73) 

The first one arises right from the definition in Eq. fl65|) and Eq. fl62l) . For the second 



formula, take /i derivatives of w^^^'^^x) in Eq. (!65|) and apply the definition of the Q 
polynomial. 



X) 
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A second property is the one which makes worthy the definition of the famihes Q^°''^\ 
for it captures the exact relation between a Romanovski polynomial and its derivative, 
and expresses it in a simple fashion. We already know that if R^n'^\x) is a Romanovski 
polynomial, so is its derivative, because it is a solution of the same hypergeometric equation 
for the parameters {a, (3 +1). So '^'^"^x ^^"-^ -^i-i'*'^'' (^) ^lust be proportional. Working 
through Eqs. ([8]), ([9]) and fl3T|) . the exact relation results: 



n(2/3 + n-l)i?i"_:^+')(x). (74) 



dx 

This equation, when expressed in terms of the Qu{x) polynomials, takes the form 

^^^^ = u{2p + u-l)Qtf\x), (75) 

which tells us that in the family Q^"''^\ each polynomial is the derivative (up to a constant 
factor) of the following polynomial in the same family. 

For a third property, let us note first that the family obeys the following differential 

recurrence relation: 

qK^(x) = s{x)^^^^^ + [2{(3 + iy-l)x + a]Qi^^^\x). (76) 

For a proof, derive the definition of Q^'^\x), Eq. fl65l) . and get 

dQi"''^\x) _ 1 dw("'^-^)(x)d'^M;(°'^)(x) ^ 1 d''+'^w^''^^\x) 



dx (w("''^-^)(x))2 dx dx'' w("'^-'')(x) dx'^+i 

The second term on the right is s~^{x)Q^"_^i\x). The first term in the right, after the 
derivation of w^°''^~'^^ (x) , gives s^^{x)[2{f3 — v — l)x + a\Q^v'^^[x). The result arises after 
reordering. Then, the substitution of the derivative formula, Eq. fl75l) . gives a three term 
recurrence relation. 

Q(";f (x) - [2(/3 + z/ - l)x + a]Q^^)(x) - ^2/3 + - 1)(1 + x^)Qt^Xx) = 0, (77) 

from which the polynomials can be efficiently generated, in contrast to the Rodrigues formu- 
las. A fourth property states that qI^'^^ (x) is the polynomial solution to the hypergeometric 
differential equation 

(1 + x2)Q"(x) + [2(/3 - u)x + a]Q'{x) + A,Q(x) = 0, (78) 
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where Xi, = — z/(2/3 — — 1). This equation is just the hypergeometric differential equation 
(1371) for the polynomial -Rt"''^ '^\^)^ which is Q^i^''^\x), so it is proved. 

The fifth property is that there exists a generating function Q^^'^\x^y) in closed form 
for the family Q.^°''^\ which is an analytic function whose Taylor expansion in the variable 
y is given by 

Q("'^)(x,y) = 5^^Q("'^)(a;). (79) 



u=0 



By substituting the definition of Q'f''^\x), Eq. ( l65i) . in previous equation and grouping 
factors and setting z = x + ys{x) we get 



1 ^ (.-x)- d- 



(80) 



which is a Taylor expansion of the function inside the derivative at the point {x + ys{x)) 
with base point x. Thus, the summation of the series is given by w^"'^) at the point B = 
X + ys{x) = X + y{l + x^). The result is 

(^ _|_ d2\/3-1 -acot-^B 

Q'"'''(-.v)= \l\J^,J,-....-^. ' (81) 



rom which numerous recursion relations, such as Eq. (1771) . may be derived as usual |17l |. 
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i8i, m. 

We now address an orthogonality property of the Qi^''^\x) polynomials. Polynomials in 
the family 2*^"'^^ with P < e — |, satisfy the following relation: 

-^Qt^Q^,^^,^ (82) 

.oo S[X)2 S{X)2 S{X)2 

where m ^ n and e = 1 if m + is odd, and e = 2 if m + n is even. This is an orthogonality 
integral between the functions Qm'^'' / s{x)'^/'^ and Q^T'^V'^l^)"^^ built on top of the poly- 
nomials in In contrast to the orthogonality relations in the 7^("'^) families, which 
are valid only for a finite subfamily of polynomials, Eq. (1821) applies to the whole family 
Q(a,/3) terms of the Romanovski polynomials in Eq. (!69l) the integral in Eq. (182|) takes 
the following form 

/^(a,/3-m)(2.) Ri^^P-rn) /^(a,/3-n) (3,) R^^^P-^) (x) —1— dx = 0, (83) 
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which can be read as orthogonahty within the infinite sequence of polynomials -R^,"'^ ''^ with 
a running parameter attached to the polynomial degree. 

In fact, Eq. (152]) stands for two different results which require separate proofs. In any 
case, since m ^ n, we can take m > n. Let us consider first the case of even m + n. Then, 
the integral of Eq. (!82|) is 

0.„,„.r!^«^«^d.. (84) 

Upon substitution of Qm'^\x) by its definition, Eq. (165|) . we get 

/CO IJTt 
s(x)^("^-'^)-ig^^)(x)^«;("'^)(x) dx. (85) 

Because m + n is even and m > n, then m — n — 2 is an even, nonnegative, integer. Thus, 
s(x)5('"^"')^^ is a polynomial of degree m — n — 2 and s(a;)5(™^")^^Q^"''^'*(x) is a polynomial of 
degree m — 2, which we call Pm-2- Then, after m—1 integrations by parts, m — 1 derivatives 
are applied to Pm-2 so it vanishes and we are left only with the boundary terms. 



0^,„ = ^(-l 



m—1 

(86) 



k=l 



d^-^Pm-2{x) d"'-''w'-'''^\x) 



For each k, the derivative of Pm-2 is a polynomial of degree m — k — 1 whereas the m — k 
derivative of is given in terms of the polynomial Q^^f^, again by its definition in 

Eq. (!65|) . Then, the k boundary term results 

d''-^Pm-2{x) d"'-''w^°'^\x) 



dx'^ 1 (l^m k 



g-«cot- ^ Xy-"'+''-'P2m-2k-l{x), (87) 



where P2m-2k-i is a polynomial of degree 2m — 2k — 1. The asymptotic behavior of this 
term at ±oo is the same as and, thus, it goes to zero if, and only if, /? < |. 

The proof of the case with (m + n) odd is similar. 



IV. RELATIONSHIP BETWEEN ROMANOVSKI POLYNOMIALS AND JA- 
COBI POLYNOMIALS 

Romanovski and Jacobi polynomials are closely related. In fact, it is common that 
Romanovski polynomials are referred to as complexified Jacobi polynomials jlO|. In 
this section we are showing which is the precise relationship between them and which is not: 
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Romanovski polynomials can indeed be obtained from a generalization of Jacobi polynomials 
to the complex plane, but not through the complexification of Jacobi polynomials, which is 
a different issue. Let us distinguish both concepts. 

For ease of reference, we recall here equations ([1]) and ([2]), which are the equations 
Romanovski polynomials and Jacobi polynomials solve, respectively. 

(1 + x^)B" + t{x)E' + \R = 0, (88a) 
{l-x'^)P" + t{x)P' + XP = 0, (88b) 

where t{x) is a polynomial of, at most, first degree. 

The argument of the complexification is based on the fact that the change from x to ix 
transforms the coefficient (1 — a;^) in Eq. fl88bl) into (l + x^), the coefficient in Eq. flHSal) . But 
caution is needed with this idea. Complexification is a transformation which takes real valued 
functions of a real variable into complex functions of a real variable. U g : U G H. H. is 
such a real valued function, defined in an open subset of the real line, we define the function 
g : W G U ^ C hj the recipe (wherever it makes sense) 

g{x)=g{tx). (89) 

The new function g may, or may not, be well defined in all the points of U and may, or may 
not, inherit the continuity and differentiability properties of g in all points of U (think, for 
an instance, of the function g{x) = (1 + x^)~^). In the case g happens to be a polynomial 
it is easy to see that both continuity and differentiability are indeed respected. Then, the 
derivatives of g and g with respect to x satisfy the following identity: 

(90) 

ne{0,l,...}. 

Complexification respects the sum and product operations, i.e., f + g = f + g and f g = fg 
(it is a ring homomorphism from a ring of real valued functions of a real variable to the 
ring of complex valued functions of a real variable). Thus, if the function g is a solution 
to a linear differential equation, then is a solution of the complexification of that linear 
differential equation. The application of this argument to the Jacobi polynomials gives the 
result that P^'^\x) = Pn"'^\ix) verifies the complexification of Eq. (I88b[) . namely 

(1 + x^){P)" + it{ix){P)' - AP = 0, (91) 
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where the prime still stands for derivative with respect to the real variable x. But equation 
(I9TII is not the same as equation fl88ap unless it{ix) is real. If we write t{x) = (3 — a — {a + 
j3 + 2)x, as is customary in the Jacobi equation (see Eq. (fT9|) ). we need (a + /? + 2) to be 
real and [jS — a) to be imaginary, which is achieved only if a and (3 are complex and (3 = a*. 
Hence we have to consider the functions P^'^\ix) with complex parameters a and (3 which 
are no longer the complexification of the classical Jacobi polynomials as described above. So, 
the complexification of Jacobi polynomials does not result in the Romanovski polynomials. 
But, even in case it did, not all the properties of Jacobi polynomials would be translated to 
properties of Romanovski polynomials: only those which made use of theorems like equation 
(!90|) . which relates the derivatives. For instance, there is no theorem relating integrals of 
a complexified function and the original function; thus, all the properties depending on 
integrations, such as the orthogonality, would have no translation to the complexified version. 

An alternative scenario is to extend the definition of Jacobi polynomials to the complex 
)lane: complex variable z, complex parameters a and (3 and, obviously, complex values. In 



20l | this definition has been successfully given as 

- i E (;; 1 1) (" t ") (1 - ^)'<i + '^2' 

or, equivalently, by the Rodrigues formula 

1 (7" 

plc^'\z) = - zr-{i + z)-^— [{I - zT+-{i + zr^] , (93) 

which are formally the same as the classical ones except now z,a,(3 G C while n G 
{0, 1, 2, . . . }. These polynomials solve the complex Jacobi ODE 

(1 - z^)P'\z) + [l3-a-{a + P + 2)z\P\z) + (a + /3 + 1 + n)nP{z) = 0, (94) 

where the prime stands now for the derivative with respect to the complex variable z. 

The specialization of the variable to the imaginary axis, z = ix, and the parameters to 
f3 = a* leaves us with Eq. fl88ap (notice the change from d/dz to d/dx gives an extra i), so 
these complex Jacobi polynomials solve the differential equation that define the Romanovski 
polynomials. One has to prove now that functions Pn"'" \z) restricted to z = ix are real 
valued or, at least, proportional to a real valued one. This is easily achieved by computing 
the complex conjugate of the function Pi"'° \ix) in terms of the definition in Eq. fl92l) 

^ ^ t (::;') (" r) (1 + ^^f<^ - (95) 
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With a change in the summation index from k to I = n — k, we get 

Pi"'"')(zx)* = (-l)"Pi"'"*)(^a;), (96) 

i.e., for even n, Pr^'" \ix) is real, while for odd n, it is imaginary. Hence, the combination 
i^-Pn"'" \ix) is a real function for all n. Finally, because the polynomial solutions of a 
hypergeometric differential equation are unique for each degree, up to a constant factor, 
we conclude that z"Pi"'° \ix) is the Romanovski polynomial of degree n with parameters 
— 253(q;) and — (3?(«) + 1). In other words, complex Jacobi polynomials do provide another 
characterization of the Romanovski polynomials via 

Rl:'^\x) = ^"Pi'-^-^"''-^+^"\^a;), (97) 

(with suitably chosen normalization constants for the Jacobi polynomials). However, this 
alternative characterization is of no help when it comes to study the orthogonality properties, 
because the orthogonality properties of the complex Jacobi polynomials are not well known. 
In the authors state some new results on the orthogonality along some particular paths 
in the complex plane. For instance, in their Eq. (4.3) an orthogonality relation is given 
along the imaginary axis, which is our case, but only for a special case demanding real, not 
integer, parameters, which is not our case. To our knowledge, at the present time there are 
no results concerning the orthogonality of these complex polynomials which would provide 
an alternative approach to the results on orthogonality stated in Subsection IIII A[ 

The argument presented here states that Romanovski polynomials are just a subset of 
complex Jacobi polynomials. Therefore it may seem that Romanovski polynomials are, 
somehow, subordinated to the Jacobi polynomials. But the whole argument could be re- 
versed if we had a definition of complex Romanovski polynomials as the one given in Ref. [2^ 
(here reproduced in Eq. (l92l) ). If that would be the case, it would not be surprising to get 
a relation of the form 

pi-'^){x) = -(z")P?"-''^'^^"+^^"''\zx), 

where, now, P^'^\x) is a real Jacobi polynomial and Pi*''" ^^"^'^'^^^''^^\ix) would be a 
complex Romanovski polynomial. But, as we do not have such a definition, this last formula 
is nothing but a conjecture. 
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V. ROMANOVSKI POLYNOMIALS IN SELECTED QUANTUM MECHANICS 
PROBLEMS 



The Romanovski polynomials are part of the exact solutions of several problems in ordi- 
nary and supersymmetric quantum mechanics. In this section we review a few prominent 
cases. The selection of the examples certainly reflects personal preferences and does not 
pretend to be complete. 

In general, the exactly soluble Schrodinger equations enjoy a special status because most 
of them describe phenomena that play a key role in physics. Suffice it to mention in that 
regard such textbook examples as the description of the spectrum of the hydrogen atom in 
terms of the Coulomb potential 1 2]| , or the description of vibrational modes in molecules and 
nuclei in terms of the Hulthen and Morse potentials 22|, 231]. More recently, exactly soluble 
potentials acquired importance within the context of supersymmetric quantum mechanics 
(SUSYQM) which considers the special class of Schrodinger equations {H{z) — E)'${z) = 0, 
with H{z) standing for the Hamiltonian (of the one-dimensional, real variable z), and E 
for the energy, which allow I24I a factorization of H{z) according to H{z) = A^{z)A~ (z) + 
Egst, and A~ (z)"^ gst{z) = 0. SUSYQM provides a powerful technique for finding the exact 
solutions of Schrodinger equations. To be specific, any excited state can be obtained by 
the successive action on the ground state, \l/gst(-2), of an appropriate number of creation 
operators, A^{z), defined in terms of the so-called superpotential, U{z), as 



Supersymmetric quantum mechanics governs a family of exactly soluble potentials (see 



Refs. |25|- [28| for details) two of which are the so-called hyperbolic Scarf and trigono- 

' n Q n 

metric Rosen- Morse potentials, that have been solved recently in [6|, [7|], [16| in terms of the 
Romanovski polynomials as discussed in the next two subsections. The third subsection is 
devoted to applications of the Romanovski polynomials in random matrix theory. 
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A. Romanovski polynomials in problems with non-central electric potentials 



The (one-dimensional) Schrodinger equation with the hyperbohc Scarf potential is 

d2 



V,{z)^[B'-A{A + l)]- 



2/i dz2 
1 



+ VJz)-E l^ z) = 



(98) 



B{2A+ l)tanh2- 



' cosh^ z ^ ' ' - j^Qg]^ ^ 
This equation appears, among others, in the problem of a particle within a non-central scalar 
potential, a result due to Ref. j^. In denoting such a potential by V{r, 6), one can make for 
it the specific choice of 

V2{e)_ 

' ' (99) 



V{r, 



Vi(r) + 



V2{e) = -hcote. 



An interesting phenomenon is the electrostatic non-central potential in which case Vi(r) is 
the Coulomb potential. The corresponding Schrodinger equation 



\ l d 1 



2/i 



d . 1 
■ sm U— + 



dr dr sin 9 89 89 sin^ 9 86'^ 



+ V{r,9) 



^{r,9,ip) 
= E^{r,9,v), (100) 



is solved in the standard way by separating variables. As long as the potential does not 
depend on the azimuthal angle, one assumes 



^(r,^,0) =7^(r)0(0)e^'^"^. 

The radial and angular differential equations for TZ{r) and O(^) are then found as 
d^n{r) 2 dn{r) 



(101) 



and 



d^Q{9) 
d9^ ^ 



+ 



r dr 



+ 



,^J<S)(9) 



n{r) = 0, 



0(9) = 



(102) 



(103) 



sm'9_ 

with /(/ + 1) being the separation constant. From now on we will focus on the second 
equation. Notice that for V2{9) = 0, and upon changing variables from 9 to cos 9, the last 
equation transforms into the associated Legendre equation and correspondingly 



0(0) ^i^p-(cos( 



(104) 
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an observation that will become important below. 

Following Ref. [9i] one begins with substituting the polar angle variable by a new variable, 
2;, introduced via 9 = f{z), with / to be determined. This leads to the new equation 



+ 



+ f'iz) cot fiz) 



d_ 

dz 



+ 



m 



'^V,ifiz)) + lil + l)--^-- 
n sm / [z) 



i){z) =0 (105) 



with f'{z) = 



dfjz) 
dz ' 



and Tpi^z) defined as ipi^z) = Q{f{z)). Next one can require that f'{z) 



approaches zero at z = like sinz, meaning, lim^^o /'(-z)/ sinz = 1, and define f{z) via 



f{z) cot f{z). 



The latter equation is solved by f{z) = 2 tan ^ e^. With this relation one finds that 

1 



sin 6' 



cos^^ 



tanh z 



(106) 



(107) 



cosh z ' 

and consequently, f'{z) = sin f{z) = sechz. Upon substituting the last relations into 
Eqs. f l99|) . and fllOSp . one arrives at 



(fij{z) 
dz^ 



+ 



1(1 + 1) 



1 2u, , 1 

5 5-0 tanh z — 

cosh z h cosh z 



m 



ip{z) = 0. 



(108) 



In taking in consideration Eqs. fl98l) . (11071) one realizes that the letter equation is precisely 
the one-dimensional Schrodinger equation with the hyperbolic Scarf potential and with 

• /(/ + 1) playing the role of -{B^ - A{A + 1))/ {n^/{2fi)), 

• m? playing the role of —En/ (/?.^/(2yu)), 

• h playing the role of —B(2A + 1). 

This equation has been solved in terms of the Romanovski polynomials in Ref. 6| upon sub- 
stituting sinh^; = x. Notice that there the weight function was defined as (l+x'^)~^e'^^'^'^ 
and the polynomials have been labeled correspondingly as Rn'''\x), following 29 1. A com- 
parison with Eq. (l25l) allows identifying p — > —(3 + 1, g — > a. In terms of the notations 



of the present work, the result of Ref. 



can be cast into the following form: 



^^[z) = Cn[l + (sinhz)2]§-4 et*'^'^ '^'°^^"i?(^°'^)(sinh 

1 



a 



-2B, (3 = -A 



E.n 



-{A-n) 



(109) 
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with Cn being a normalization constant. Back to the 6 variable and in making use of the 

def 

equality x := sinh z = — cot 6, we find 

e(^) = ^„(sinh^(-cot^)) = Cn[l + (cot ^)'] 2-3 et*"'^"'(-^°*')i?^'^)(-cot^), (110) 

showing that the angular part of the exact solution to the non-central potential under con- 
sideration is defined by the Romanovski polynomials. In turning off the non-central piece 
of the potential, the angular part of the solutions will become the standard spherical har- 
monics, Yi^{6, (p) = Pl^{0) e*™''^, which will produce a relationship between the Romanovski 
polynomials and the associated Legendre functions, an issue to be considered in more detail 
at the end of this section. 

Now, in accord with the theorem on the finite orthogonality of the Romanovski polyno- 
mials in Eq. (!55|) . also only a finite number of eigen-wave functions to the hyperbolic Scarf 
potential appears orthogonal, 

Mz)i^m{z)dz = CmCn \ W^''^<'\x)Rt^<'\x)R^:^<'\x) = Smn, 

Joe (111) 

m + n < 2A . 

This finite orthogonality refiects the finite number of bound states within the potential under 
consideration. 

Next, it is quite instructive to consider the case of a vanishing ¥2(6), i.e. 6 = 0, and 
compare Eq. (11081) to Eq. fl98|) for = 0. From now on we will give all quantities in units 
of h = 1 = 2fi. In this case 

• Eq. fllOSp reduces to the equation for the associated Legendre polynomials, P;'"(cos6'), 

• / becomes A, 



rn? becomes (/ — n)^. 



• Eq. dHBD produces ^^'^\x) as part of its solutions, 

which allows one to relate n to / and m as m = I — n. In taking into account Eq. (I104p 
together with cot 6 = ~ sinh z provides the following relationship between the associated 
Legendre functions and the Romanovski polynomials 

pricose) = const[i + {cote)']"^R'~^f\-cote) , 

(112) 

m + / = n e {0, 1, . . . ,/}. 
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(114) 



In substituting the latter equation into the orthogonahty integral between the associated 
Legendre functions, 

"1 

pr (cos e)PiT {cos 6) d cos e = 0, 

(113) 
we find 

Jjl + (cot^)^)-^<'i-'Vcot^)<t'Vcot^)dcos^ = 0, 

I ^ I'. 

The latter relationship amounts to the following orthogonality integral 
/""(I + (sinh^)2)-^i?i°'^"'\sinhz)i?j;!'^"''\sinh2)(sech z)Mz = 

J —oo 

r y/ti;(o.|-0 (x)i?i°'^~'^ (x) (x)i?i°'^~''^ (x) dx = 0, 

x = sinh^, l-n = l' -n' = m>0, l^l'. (115) 

Careful inspection shows that this equation is nothing but a particular case of the orthogo- 
nality relation in the family of polynomials Q^°''^\ established in Eq. ( l82l) and translated to 
the Rt''^^ notation in Eq. (1831). 

A further example is given by the Klein-Gordon equation with equal scalar and vector 
potentials. It has been shown in Ref. that the former can be reduced to the corresponding 
Schrodinger equation. Therefore, in case one uses the hyperbolic Scarf potential in the above 
Klein-Gordon equation, one will face again the Romanovski polynomials as part of its exact 
solutions. 



B. Romanovski polynomials in quark physics. 

The interaction of quarks, the fundamental constituents of the baryons, are governed by 
Quantum Chromodynamics (QCD) which is a non-Abelian gauge theory with gauge bosons 
being the so called gluons. QCD predicts that the quark interactions run from one- to many 
gluon exchanges over gluon self-interactions, the latter being responsible for the so-called 
quark confinement, where highly energetic quarks remain trapped but behave as (asymptot- 
ically) free particles at high energies and momenta. The QCD equations are nonlinear and 
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FIG. 1: The non-central potential V{r,6), here displayed in its intersection with the x = plane, 
i.e. for r = y^y^ + 2;^, and 6 = tan^^ ^. The polar angle part of its exact solutions is expressed in 
terms of the Roman ovski polynomials. 



complicated due to the gluonic self-interaction processes and their solution requires employ- 
ment of highly sophisticated techniques such as discretization of space time, so-called lattice 
QCD. Lattice QCD calculations of the properties of hadrons, which are all strongly interact- 
ing composite particles, predict a linear confinement potential with increasing energy. The 
one-gluon exchange potential, which is Coulomb-like, ~ 1/r, adds to the linear confinement 
potential, ~ r, a combination that is believed to provide the basic properties of two-body 
(mainly quark-antiquark) systems as concluded from quark model calculations. In contrast 
to this, three quark systems (baryons) have been believed for a long time to involve more 
complicated interactions depending on the internal quark degrees of freedom such as their 
spins, isospins, flavors, and combinations of them, while the potential in coordinate space 
has been considered of lesser relevance and modeled preferably by means of the harmonic 
oscillator. However, in so doing, one encounters the problem of a serious excess of predicted 
baryon excitations in comparison with data 30|] (so called "missing resonances"). 

More recently, the structure of the baryon spectra has been re-analyzed in Refs. 3l| with 
the emphasis on the light quark resonances. The result was the observation of a striking 
grouping of resonances with different spins and parities in narrow mass bands separated 
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by significant spacings. More specifically, it was found that to a very good accuracy, the 
nucleon excitation levels carry the same degeneracies as the levels of the electron with spin 
in the hydrogen atom, though the splittings of the former are quite different from those of 
the latter. Namely, compared to the hydrogen atom, the baryon level splittings contain, in 
addition to the Balmer term, also its inverse but of opposite sign. The same was found to 
be valid for the excitation spectrum of the so called A(1232) particle, the most important 
baryon excitation after the nucleon. The appeal of these results lies in the fact that no 
state drops out of the systematics, on the one side, and that the number of "missing" 
states predicted by it is significantly less than within all preceding schemes. The observed 



degeneracies in the spectra of the light quark baryons have been attributed in Ref. [32(] to 
the dominance of a quark-antiquark configuration in baryon structure. Within the light of 
these findin gs, the form of the potential in configuration space acquires importance anew. 
In Refs. the case was made that the trigonometric Rosen- Morse potential provides 

precisely degeneracies and level splittings as required by the light quark baryon spectra. The 
trigonometric Rosen- Morse potential in the parametrization of Ref. [16| reads: 

VtRM{z) = -2bcotz + l{l + l)-^, (116) 

sm z 

with / standing for the relative angular momentum between the quark and the di-quark in 
units, as usual, of h = 1 = 2/i, and 2; = ^ is a dimensionless variable built with a suited 
length scale d. 

The reason for the success of this potential in quark physics is that it captures the essential 
traits of the QCD quark-gluon dynamics in interpolating between the Coulomb potential 
(associated with the one-gluon exchange) and the infinite wall potential (associated with the 
trapped but asymptotically free quarks) while passing through a linear confinement region 
(as predicted by lattice QCD) (see Fig. 2). It is quite instructive to perform the Taylor 
expansion of the potential of interest, 

viz)tRM ^ -7 + y ^ + + ^T5 + - ^^^^^ 

This expansion clearly reveals the proximity of the cot term to the Coulomb-plus-linear 
confinement potential, and the proximity of the csc^ term to the standard centrifugal bar- 
rier. The great advantage of the trigonometric Rosen-Morse potential over the linear-plus- 
Coulomb potential is that while the latter is neither especially symmetric, nor exactly soluble, 
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FIG. 2: The trigonometric Rosen- Morse potential (solid line) and its proximity to the Coulomb- 
plus-linear potential as predicted by lattice QCD (dashed line) for the toy values / = 1,6 = 50 of 
the parameters. 



the former is both, it has the dynamical 0(4) symmetry (as the hydrogen atom) and is ex- 
actly soluble. The exact solutions of the, now three dimensional, Schrodinger equation with 
VtRAiiz) from Eq. f lll6p have been constructed in on the basis of the one-dimensional 



solutions found in 



and read: 



Mcot-'x) = (l + a;2)-^e-™^°*"'^-C^^"^'^+'™\x), (118) 

with X = cotz. The C polynomials from j?] are Romanovski polynomials but with parame- 
ters that depend on the degree of the polynomial. The following identification is valid: 

d"^"^'^'"''™^(x)^4-''^")(x), 

an= , , , , , /3„ = -(n + O + 1, n G {0,1,2,...} 
n + I + 1 

The Schrodinger wave functions are orthogonal, as they are eigenf unctions of a Hamilton 
operator. Their orthogonality integral (here in z space) reads 

dz i)n{z)i)n'{z) = 5nn' ■ (120) 







The orthogonality of the wave functions ipniz) implies in x space orthogonality of the 
rI^"'^"\x) polynomials with respect to w'--°'"'^"\x)^ due to the variable change. As long as 
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dx 

""^ dx 



dcot ^x _ -1/(^1 _|_ a;2-j _ then the orthogonahty integral takes the form 

To recapitulate, the Romanovski polynomials have been shown to be important ingredients 
of the wave functions of quarks in accord with QCD quark-gluon dynamics. 

C. Romanovski polynomials in random matrix theory 



Random matrix theory was pioneered by Wigner [33[] for the sake of modeling spectra 
of heavy nuclei which are characterized by complicated interactions between large numbers 
of protons and neutrons. Wigner's idea was to limit the infinite dimensional Hamiltonian 
matrix in configuration space to a finite, real quadratic (A^ x A^), and symmetric, matrix 
with elements being chosen at random from a suitable probability density distribution, say, 
the Gaussian one. Along this line one can then model the densities of the nuclear states 
as averages over the weighted sets of matrices. The advantage of this method is that as 
A^ — > 00, the (normalized) eigenvalues of any randomly chosen matrix approach the limits 
of the corresponding system averages, much like the general limit theorem. The probability 
density distribution (p.d.f.) of the eigenvalues of the Gaussian ensemble of random matrices 



is given by (the presentation in this section closely follows Ref. |34|): 



^e-^^U^ n l^^-^.l' (122) 

where Cn is the normalization constant and Aj are the eigenvalues. Besides the Gaussian 
random ensemble there are other random matrix ensembles under consideration in quantum 
physics such as the circular Jacobi ensemble and the Cauchy ensemble, and precisely these 
are of interest in the present section. However, their definitions require us to go beyond 
the ensembles of random real matrices and consider matrices with complex entries. To be 
more specific, one considers random ensembles composed by symmetric, unitary matrices, in 
which case the theory is not developed from an explicit distribution density function for their 
elements but rather from the requirement of the existence of a certain appropriate uniform 
measure. The random unitary matrix ensemble is special not only because it forms a group 
but mainly because this group is compact and allows for the definition of the so called Haar 
volume, which then provides the uniform measure on the space as required above. 
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Definition: The circular unitary ensemble is the group of unitary matrices U 
endowed with the volume form (rf//U) = ^ (U^U) = id'M.2 with M2 hermitian. 

The eigenvalues of the circular ensembles are confined to the unit circle, i.e. to = e*^^ 
with — TT < 9j < n. The associated probability density distribution of the eigenvalues is then 
given by 



^ n 



C 

l<j<k<N (123) 
-n < 9i < n . 

More generally, an ensemble of unitary and symmetric matrices has an eigenvalues p.d.f of 
the form (in the notations of Ref. ^34] ) 

N 



n«^2(^o n 



\Zk — z 



'=1 ^<j<k<N (124) 

where W2{zi) is a specific weight function. The circular Jacobi ensemble is specified by 

^2(2) = |1 - zp". (125) 

A relevant research goal in quantum physics is finding the spacings in the spectra of the 
circular Jacobi ensemble. Compared to the state densities, the calculation of gap proba- 
bilities in the spectra deserves special efforts. In this section we review briefly the concept 
for the calculation of gap probabilities in the circular Jacobi ensemble by means of the so 
called Cauchy random matrix ensemble, a venue that will conduct us one more time to the 
Romanovski polynomials . 

To begin with, one considers the mapping 

= . (126) 

which maps each point A on the real line to a point 9 on the unit circle (measured anticlock- 
wise from the origin) via a stereographic projection. Changing correspondingly variables in 
Eqs. f ll24l) - fll25p amounts to the following eigenvalue p.d.f.: 

N 

n(i+Ai)-^-°' n i^^-^^-i' 

1=1 l<j<k<N (127) 

Xj E ( — 00, +00) . 
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As long as one recognizes in the weight function the Cauchy weight, the random unitary 
matrix ensemble generated this way is termed the Cauchy ensemble. On the other hand, 
a comparison with the weight function of the Romanovski polynomials reveals the Cauchy 
weight as w^^'~'^~°''^^\x), an observation that will acquire a profound importance in the 
following. 

Back to the main goal, the gap probability, or better, the probabihty for no eigenvalues 



in a region J, denoted by -£^(0, /), and for the case of any ensemble is now calculated as 34 1 



°° (-1)" f f I 

E{0,I) = 1 + ^ j— / dxi... / dxndet'^ \/ W2{xi)pi{xi)^J W2{xj)pi{xj) , (128) 

n=l "'^ 1=0 

where pi{x) with / = 0,1,2,... stand for the orthogonal polynomials associated with the 
weight function ^2(2;). In other words, knowing the orthogonal polynomials is crucial for 
the calculation of gap probabilities in any random matrix ensemble. In the specific case under 
consideration, one seems to have two options in the choice for those polynomials, Romanovski 
versus Jacobi polynomials in accord with their relationship established in Eq. (!97|) . The 
choice is clearly in favor of the Romanovski polynomials because the formalism developed 
for calculating £"(0, /) (see Ref. jssi for details) is based on real coupled differential equations. 
In choosing the Romanovski polynomials one, to speak with the authors of Ref. |3^] , avoids 
the clumsy and unnecessary work of recasting the formalism on the circle, i.e. in terms of 
the complexified Jacobi polynomials. In summary, the Romanovski polynomials (termed 



Cauchy weight polynomials in Ref. 



34| ) provide a natural and comfortable tool for finding 



all the results for the circular Jacobi ensemble from those of the Cauchy ensemble. 

VI. ON THE ORTHOGONALITY RELATIONS 

We have shown in the previous sections not one but several different orthogonality rela- 
tions among the Romanovski polynomials. In this section we comment on this issue. 

First we have shown, in Eq. 0551) . a finite orthogonality in the family 7^("'^) (that of 
polynomials with fixed parameters a and (3). It is the equivalent relation to the well known 
orthogonality of the Hermite, Laguerre and Jacobi polynomials, except that in these cases, 
there is complete orthogonality. The finite orthogonality, however, is required as such in 
the solution of the wave eigenfunctions of the hyperbolic Scarf potential, see Eq. (11 111) in 



subsection IV Al or Ref 



•El: 



in this case only a finite number of states are bounded; precisely 
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those which are normahzable. The orthogonahty relation in the family, Eq. (I82l) . 

is a complete, not finite, orthogonality: it is valid for all the polynomials in the family. 
The difference with the previous is that the Q^"'^) family is made up with Romanovski 
polynomials with one fixed parameter but the other running attached to the degree. This 
different orthogonality also find its application in, for instance, Eq. (11151) in subsection IV A[ 
Finally another physics problem, the eigenf unctions of the trigonometric Rosen-Morse 
potential studied in [7] and jlO] and revised in Subsection IV Bt has given rise to yet another 
orthogonality relation: the one in Eq. (11211) . which is very similar to the orthogonality 
in the family, but not equivalent. The polynomials involved in Eq. (I12ip have both 

parameters, a and /3 running with the degree, as shown in Eq. (I119p . This last orthogonality 
is proved not directly as the others, but by means of the Schrodinger equation where it comes 
from: as the functions involved are the eigenfunctions of a self-adjoint operator, they are 
orthogonal. Here, thus, it seems as if the Schrodinger equation carefully chooses, from 
the set of all Romanovski polynomials, another family with a special combination between 
parameters and degrees such that another orthogonality relation surfaces. This kind of fine 
tuned combination of parameters is not completely new. Here is a well known instance: the 
radial part of the well-known solution of the hydrogen atom, which is given by 



nni{xn) = Nni^^--K'^iU{^n). A = 2/ + 1, X, = a„r, (129) 



where x„ is the dimensionless but n dependent variable (see also Problem 13.2.11 in Ref. [19|), 
while r is the radial one. Here L^m^^ [x) is the generalized Laguerre polynomial of degree 
m as introduced after Eq. f|T7|) . Notice that Laguerre polynomials of different degrees in 
Eq. (11291) emerge within different potential strengths, Ze^/a„, and, henceforth, the orthog- 
onality relation given by the Schrodinger equation 

oo - - 
V^^n y/Xn' ^ ' (130) 

f3k = 2k + 1, m{n,k) =n - k ~1, n^n', Xn' = —Xn, 

is not equivalent to the orthogonality given by the weight function, Eq. ( |T8|) . here restated 
for a = 1 

x^ ll)^^^ {x)x^ L^lf\x)dx = 0, 

(131) 

m 7^ m', j3 > 0. 
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In particular, notice that Eq. (11311) . when /5 G N, is recovered from Eq. (11301) in the case 
/' = /, but for /3 ^ N both formulas are completely different. 

In the Introduction we said that, perhaps, the lack of general orthogonality of Romanovski 
polynomials has been seen as a weakness and because of it they have not attracted as much 
attention as the classical orthogonal polynomials. Now we have shown that, far from being 
a weakness, the various orthogonality relations of Romanovski polynomials give them new 
appealing properties which widens their possible applications. 

VII. CONCLUSIONS 

We have presented a fairly complete description of the Romanovski polynomials as so- 
lutions of the hypergeometric differential equation ([T]), properties derived from it and some 
applications, with the following prominent items: 

1. We have described a complete classification of the hypergeometric differential equations 
in order to place Eq. ([1]) in its proper context. 

2. We have described completely, in Eq. fl28|) the polynomial solutions to Eq. ([1]), which 
are the Romanovski polynomials. We have also stated some known and some new 
properties of these polynomials. We have proposed different partitions of the set of 
all Romanovski polynomials into families which allows one to express the plethora of 
properties in a simpler and more ordered form. This approach can be applied as well to 
the other four classes of polynomial solutions of hypergeometric equations: Hermite, 
Laguerre, Jacobi and Bessel. 

3. In particular we have stated exact results about several orthogonality relations among 
the Romanovski polynomials. We have shown that a family of Romanovski polynomi- 
als, solutions to the same hypergeometric equation, is not completely orthogonal, but 
exhibits a finite orthogonality, Eq. ( 155|) . However, we have found two other orthog- 
onality relations, in families with running parameters (attached to the degree of the 
polynomial) which provide infinite orthogonality, Eqs. (182|) and (11211) . 

4. The relationship between Romanovski polynomials and Jacobi polynomials has been 
precisely stated: Romanovski polynomials cannot be obtained as just a complexifi- 
cation of Jacobi polynomials (i.e., change x by ix), but they can be realized as a 
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particularization of complex Jacobi polynomials (an extension to the complex plane 
with complex parameters). Yet, despite this relation, these complex Jacobi polyno- 
mials are not completely understood so, for instance, the orthogonality properties of 
Romanovski polynomials cannot be derived, at the present time, from properties of 
the Jacobi polynomials. 

5. We have presented three instances of the use of Romanovski polynomials in actual 
physics problems. In particular the polynomials introduced in [?] and 6] are recognized 
as Romanovski polynomials. The orthonormality relations shown in these references 
are explained in this context. 
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